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Abstract 

We calculate the two-pion wave function in the ground state of the 1 = 2 5-wave system and find 
the interaction range between two pions, which allows us to examine the validity of the necessary 
condition for the finite-volume method for the scattering length proposed by Liischer. We work in 
the quenched approximation employing a renormalization group improved gauge action for gluons 
and an improved Wilson action for quarks at 1/a = 1.207(12) GeV on 16^ x 80, 20^ x 80 and 
24^ X 80 lattices. We conclude that the necessary condition is satisfied within the statistical errors 
for the lattice sizes L > 24 (3.92 fm) when the quark mass is in the range that corresponds to 
mf = 0.273 — 0.736 GeV^. We obtain the scattering length with a smaller statistical error from 
the wave function than from the two-pion time correlator. 
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I. INTRODUCTION 


Calculations of the scattering length and the phase shift represent an important step for 
expanding our understanding of the strong interaction based on lattice QCD to dynamical 
aspects of hadrons. For the simplest case of the two-pion system in the 1 = 2 S'-wave system, 
the scattering length has been ealeulated in MBiBBllEa and the pioneering study 
of the phase shift was made by Fiebig et al. |l0| using the two-pion effective potential. We 


irect calculation o 
and full QCD 


11 


the phase shift without recourse to the effective potential 


12j |. Kim rep orted on preliminary results of the phase shift 


presented a c 
in quenched 

on G- and Ff-periodic boundary lattices 

The calculation of the scattering length and the phase shift usually employs the hnite- 
volume method of Liischer, in which the scattering phase shift is related to the energy 


15, 


16, 


17| |. In previous application of the formula the 


eigenvalue on a hnite volume 14, 
energy eigenvalues were calculated from the asymptotic time behavior of the two-pion time 
correlator. 

The derivation of Liischer’s formula assumes the condition R < L/2 for the two-pion 
interaction range R and the lattice size L, so that the boundary condition does not distort the 
shape of the two-pion interaction. In the studies to date, however, calculations were carried 
out without verifying this necessary condition, but simply employing a few lattices having 
different sizes and extrapolating to the inhnite volume limit assuming simple functions, 
such as an inverse power of the lattice extent, for finite volume corrections. In the absence 
of theoretical justification, however, such assumptions would cause ambiguities and it is 
important to examine the validity of the condition in lattice simulations for reliable results. 

In this work we restrict ourselves to the ground state of the 1 = 2 S'-wave two-pion 
system and calculate the two-pion wave function. We investigate the two-pion interaction 
range and the validity of the necessary condition for Liischer’s formula. We attempt to 
extract the scattering length directly from the wave function and compare it with the more 
conventional result from the two-pion time correlator. 


We refer to the work for two-dimensional statistical models in 


For the two-pion 


system of QCD a similar idea was discussed in 0. We also quote Yamazaki who pre¬ 
sented preliminary results for the wave function for the four-dimensional Ising model and 
the scattering phase shift therefrom j0. 

This paper is organized as follows. In Sec. O we give a brief review of the derivation of 


Liischer’s formula jl2l with emphasis on the role of the condition R < L/2. The calculational 
method of the wave function and the simulation parameters are given in Sec. lIIIl In Sec. IIV Al 
we present the wave function and estimate the two-pion interaction range for a 24/ lattice. 
In Sec. lIVRI we calculate the scattering length from the wave function and compare it with 
that from the two-pion time correlator. Our investigations on 20^ and 16^ lattices are given 
in Sec. IIV Cl where finite volume effects on the scattering length are examined by comparing 
three lattice volumes. Our conclusions are given in Sec. 0 Preliminary reports of the present 
work were presented in 
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II. LUSCHER’S FORMULA 


We briefly review the derivation of Liischer’s formula, with emphasis on the role of the 
condition for the two-pion interaction range. The formula [ij, [13 was rederived using an 


effective Schrddinger equation for the two-dimensional scalar hied theory in 1^, and for 
the four-dimensional case in which is discussed here. Another approach based on the 
Bethe-Salpeter wave function in quantum held theory [l^ is discussed in Appendix. 

The static two-pion wave function k) with the energy eigenvalue E = 
in the center of mass system on a hnite periodic box of volume satishes the ehective 


Schrddinger equation [l4, 


16 


(A -h k) = J d^y Uk{x, y)(j){y; k) , 


( 1 ) 


where x and y are the relative coordinate of the two pions. Uk{x, y) is the Fourier transform 
of the modihed Bethe-Salpeter kernel for the two-pion interaction on the hnite volume 


14 


and is related to the oh-shell two-pion scattering amplitude (see AppendixEj). It is generally 
non-local and depends on the two-pion energy. It should be noticed that in m is not a 
square of a 3-dimensional momentum vector but dehned from the energy by k"^ = A^/4 —m^. 
It may take a negative value in some cases. We call k “momentum” for simplicity in this 
paper, however. 

In the derivation of Liischer’s formula it is assumed that the two-pion interaction range 
is smaller than one half the lattice extent, he. there exists the distance R < T/2, where the 
wave function satishes the Helmholtz equation : 


(A -I- k) = d for X e Vr , 


( 2 ) 


with 

Vr = { X \ \x + nL\ > i? , n G } . (3) 

Next we consider solutions of 0. In general k‘^L‘^/{2nY can be integer or non-integer. 
The former case is called “singular-value solutions” in (see also 0]). The appearance 
of these solutions is not generic. It occurs only for some specihc lattice volumes or particu¬ 
lar cases of the two-pion interaction. For the two-pion ground state there is an important 
singular-value solution that has = 0, which, however, exists only on specihc lattice vol¬ 
umes or for the vanishing scattering length. If this solution appears in numerical simulations, 
the two-pion time correlator should behave as 


(0|7r(t)7r(f)7r(0)7r(0)|0) 

(0k(f)vr(0)|0)2 


const. , 


( 4 ) 


for a large t. However, such a time behavior has not been observed in previous studies or in 
our numerical simulations. Thus such a case hardly occurs for the ground state. The formula 
for the singular-value solutions was derived in Ref. ^3) but we consider only non-integer 
value solutions in this paper. 
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General solutions of the Helmholtz equation m can be written by 


0 (f; k) = vim{k) ■ Gi„i{x] k) , 

/=0 m=—l 

with vimik). Gim{x] k) is given from the periodic Green function : 

.271 


G{x;k) = 2 ^ 1,2 


^ip-x 


r = { p \ p = n— , n G Z'^} 


as 


( 5 ) 


( 6 ) 


( 7 ) 


Gim{x;k) = v^yi^{V)G{x;k) , 

where yim{x) is a polynomial related to the spherical harmonics through yim{x) = x^-Yimift^) 
with Ga; the spherical coordinate for x and x = \x\. The convention of Y;m(G) is that of 22| . 
as is adopted in [l^. It then follows that Goo{x] k) = G(x; k). 

In general we can expand the solution of the Helmholtz equation in terms of the spherical 
Bessel ji{kx) and Neumann functions ni{kx) for R < x < L/2 as 


(j){Y,k) = Y Y y^yimikl^)(ai{k) ■ ji{kx) + Pi{k) ■ ni{kx) 

1=0 m=—l 

where the conventions of ji{x) and ni{x) agree with those in 


( 8 ) 


22 


and [l^ . 


The expan¬ 
sion coefficients ai{k) and f3i{k) yield the scattering phase shift in the inhnite volume as 
tan Si{k) = j3i{k)/ai{k). In particular, for the ground state the momentum k"^ is very small 
and the S'-wave scattering length oq is given by tan5o(^) = Po{k)/ao{k) = aok + O(fc^). 

For the wave function (jH)), Pi{k) and ai{k) are geometrically related, because they can be 
expressed in terms of the expansion coefficients for Gim{x;k). The expansion of G{x;k) is 
given by 


oo m=l 

G{x]k) = —no{kx) + Y Y '/^yim{^x)gim{k)ji{kx) , 

Z=0 m=-l 


where 


L^^rP ^ 


277 

T = { p \ p = n— , n E Z^} 
Jj 


( 9 ) 


( 10 ) 


with spherical coordinate Vtp for p. (The function gim{k) differs from gim in (3.31) of Ref. [l^ 
only by the normalization as gim{k) = l/\/4^ ■ gim)- The explicit expansion for Gim,{x; k) 
with general I and m is not needed. Note that the indices I and m are not labels of the 
angular momentum, nor Gim{x] k) is the eigenstate of the angular momentum labeled by I 
and m. Actually it includes Y)/m'(Gx) with P ^ I and m' 7 ^ m. Also note that Gim,{x]k) 
contains ji/{kx) for a range of I' and nii{kx) with only V = 1. These are easily known from 
(|7|) and (jni). 

In this work we consider only wave functions in the representation of the cubic group, 
which equals S'-wave up to angular momenta I > 4. The wave function for x eVr can be 
expressed as 

0(T;/c) ='i;oo(^) ■ G'(T;/c)-F n4o(/c) ■ ^X!<^4 o(R-[^]; H-, (H) 


48 


R 
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where a vector operation R represents an element of the cnbic gronp which has 48 elements. 
The terms with I > 6 are neglected, and the I = 4 terms with m 7 ^ 0 do not appear since 
they either vanish (for \m\ 7 ^ 4) or reduce to G 4 o(T; k) (for \m\ = 4). If the scattering phase 
shift 6i{k) with Z > 4 is negligible in the energy range under consideration, (5i{k) ~ 0 for 
Z > 4 in (jHI). This means that v^o{k) ~ 0 in m, and thus 

(j){x;k) voo{k) ■ G{x;k) , (12) 

because Gim{x]k) contains niikx). This expectation is supported for the ground state by 
our numerical simulations. 

Finally we obtain Liischer’s formula between the F-wave scattering phase shift and the 
allowed value of Zc^ by comparing the coefficients of jo{kx) and no{kx) in 0 , 

1 ao{k) goo{k) 47r 1 ^ 1 

tan5o(Zi;) po{k) Zc/(47r) k 

where the region T is dehned by ca and oa is assumed. 

We remark that contaminations from inelastic scattering are likely negligible for the 
ground state, although they may become signihcant for momentum excitation states, whose 
energies are close to or above the threshold of inelastic scattering, E ~ 


III. METHOD OF CALCULATIONS 

A. Calculation of the wave function 

Our dehnition of the two-pion wave function is 

(j){x-,k) = (0| Q{x,t) | 7 r’'“ 7 r’'“; Zc) . , (14) 

where | 7 r+ 7 r+; k) is the ground state with energy 2Ek = A factor exp{2Ekt) is 

introduced to compensate the t dependence. The operator fl{x,t) is an interpolating held 
for the two-pion, which is dehned by 



where 7r+(x, f) = u{x,t)'y^d{x,t) is an interpolating operator for 7r+ at {x,t) and a vector 
operation R represents an element of the cubic group which has 48 elements. The summation 
over R and X projects out the Aj^ sector of the cubic group and the zero total momentum 
state. The wave function in m is the Bethe-Salpeter wave function projected to the A^ 
sector. It has the same properties as those of the wave function in Sec. E] and we can 
derive Liischer’s formula from it. Details are discussed in Appendix ^ The wave functions 
at all positions x are not independent. The number of independent position vectors is 
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N{N + 1)(A^ + 2)/6 with N = L/2 + 1 owing to the periodic boundary condition 0(x; k) = 
0(x + nL; k) with n ^7? and the invariance under the cubic group 0(5;; k) = 0(R[5f]; k). 

In order to calculate the wave function we construct the correlator : 


+ 1) = (0| n{x,t) W{to) W{to + 1) |0) , 


( 16 ) 


where W(t) is a wall source at time t dehned by 


H H d(x, t)'iiiu(g, t) , 


(17) 


which is used on conhgurations hxed to the Coulomb gauge. The two wall sources are placed 
at different time slices to and to +1 to avoid contaminations from Fierz-rearranged terms j^. 
Neglecting contributions from the momentum excitation states in the large t region, we 
obtain the wave function at the time slice t as 


0 (f; k) = 


F^^{x,t]to,to + 1 ) 

(xo, i; to,h + 1) ’ 


( 18 ) 


up to the overall constant, where xo is the reference position. 

We try to extract the energy eigenvalue of the two-pion and the scattering length from 
the wave function 0(x;/i;), but, for comparison, we also estimate them from the two-pion 
time correlator : 


P'-irnit) toi to T 1) ^2 ^ ( T07r(x, f, foi ^0 T 1) . 

X 

The single pion time correlator computed with the aid of the wall source, 

K{t;to) = ^^(0|7r+(x,f)hF(fo)|0) , 

X 

is used to construct the normalized two-pion correlator : 

-f^TTTT (t) tO: ^0 d“ 1) 


( 19 ) 


( 20 ) 


fl(() = 


F-nit') ^o)-Cr(C to + 1) 


( 21 ) 


In the absence of the singular-value solution that belongs to = 0, this behaves for a large 
t as 

R(t) = C ■ ^ ^22) 


where C is a constant and 


AVF = 2Ek — 2m^ = 2,^ ml. + kt^ — 2m„ 


(23) 


is the energy shift due to the two-pion interaction on the hnite volume. The momentum k"^ 
is calculated from AhF, and it can be used to estimate the scattering length via Liischer’s 
formula. 
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B. Simulation parameters 


Our simulation is carried out in quenched lattice QCD employing a renormalization group 
improved gauge action for gluons, 


^G = 


/3 



(24) 


The coefficient Oi = 
group analysis 


-0.331 of the 1x2 Wilson loop is hxed by a renormalization 


23| . and Oq = 1 — SOi = 3.648 of the 1x1 Wilson loop W 


lXl| 


X, 


by the 

normalization condition, which dehnes the bare coupling (3 = Q/g"^. Our calculation is carried 
out at /9 = 2.334. Gluon conhgurations are generated with the 5-hit heat-bath algorithm 
and the over-relaxation algorithm mixed in the ratio of 1 : 4. The combination is called a 
sweep and physical quantities are measured every 200 sweeps. 

We use an improved Wilson action for the quarks 2^ with the clover coefficient Csw 
being the mean-held improved choice dehned by 


Gsm/ = = (1 - 0.8412//3)-^/^ = 1.398 , (25) 

where is the value in one-loop perturbation theory [^. Quark propagators are solved 

with the Dirichlet boundary condition imposed in the time direction for gauge conhgurations 
hxed to the Coulomb gauge. The wall source dehned by m is set at to = 12 which is 
sufficient to avoid ehects from the temporal boundary. 

The lattice cutoh was estimated as 1/a = 1.207(12) GeV (a = 0.1632(16) fm) from 
rUp j^. The lattice sizes (the numbers of conhgurations in parentheses) are 16^ x 80 (1200), 
20^ X 80 (1000) and 24^ x 80 (506), which correspond to the lattice extent 2.61, 3.26 and 
3.92 fm in physical units. Five quark masses are chosen to give = 0.273, 0.351, 0.444, 
0.588 and 0.736 GeV^. The numbers of positions that give independent wave functions are 
165, 286 and 455 for L = 16, 20 and 24. 


IV. RESULTS 

A. Wave functions 

The two-pion wave function calculated on the 24^ lattice is exemplihed in Fig. ^ on the 
{t,z) = (52,0) plane for = 0.273 GeV^, with the reference position in (fT^ hxed at 
To = (7,5,2) {xq = |To| = 8.832). The statistical errors are negligible in the scale of the 
hgure. 

The same wave function is shown in Fig. |21as a function of x = |x| for independent data 
points. The branching of the curve seen in the hgure indicates that the wave function does 
not represent a pure S'-wave. This can be understood by the consideration in what follows. 
Let f{x) be a function depending only on x = |x| for x = [—L/2, L/2]^. The hrst derivative 


7 









is given by 


( 26 ) 


■ 

Thus /(x) satisfies the boundary condition only when d/(x)/dx = 0 at the boundary where 
at least one component of the vector x takes ±L/2. This also means d/(x)/dx = 0 for 
X > L/2 from symmetry under the cubic group. The wave function for the scattering 
system generally does not satisfy this. Hence it cannot be a pure S'-wave function, but 
receives contributions from the states with angular momenta / > 0. We expect that the 
wave function that belongs to the representation contains ji{kx) with I > 4, but not 
ni{kx) with I > 4, because 6i{k) is small for I > 4 for the two-pion ground state. This is 
supported by our results as shown latter. 

We now consider the two-pion interaction from the ratio : 


V(x; k) = 


A(j){x; k) 


0 (f; k) 

Here we adopt the naive numerical Laplacian on the lattice, 


(27) 


= IZI/(2 + A) - S(x - A) + 2/(x)l , (28) 


since k‘^ is very small and the choice of the numerical Laplacian is not important for a 
large x. Away from the two-pion interaction range, i.e. x > R, we expect that V{x;k) is 
independent of x and equals to —k"^. In Fig. iniH(x; k) is plotted for the same parameters 
as for Fig. ^ The statistical errors are again negligible. We find a very clear signal and 
H(x; k) seems to be constant for x > 8. We observe a strong repulsive interaction at the 
origin consistent with the negative scattering length of the 1 = 2 two-pion system. 

The time dependence is shown in Fig. E]for the same parameters, where the abscissa is 
X = \x\ and only independent data are plotted. We draw a line of —estimated from AW 
of ()22|1 using the normalized two-pion time correlator R{t). We see that H(x; k) approaches 
a value consistent with —for a large x. V (x; k) depends only on x and does not depend 
on t in so far as f > 44. We also confirm that the wave function 0(x; k) does not vary when 
t is varied for t > 44. 

We now consider the two-pion interaction range R. In quantum field theory the wave 
function does not strictly satisfy the Helmholtz equation 0 even for the large x region in 
a large volume lattice. Hence, with k‘^ obtained from the asymptotic time behavior of the 
two-pion time correlator, V (x; k) + shows a small tail at a large x. 

We may take the wave function as satisfying the Helmholtz equation, if V (x; k) + k"^ is 
sufficiently small compared with k"^. In the present work we take an operational definition 
for the range R as the scale, where 


_ _ l/(x; k) + _ (A + fc^)0(x; k) 

’ k^ k‘^(j){x-, k) 


(29) 
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gets small enough so that it is buried into the statistical error, with the expectation that the 
systematic error for the scattering length from the interaction of tails of the wave function 
becomes smaller than statistical errors in the resulting scattering length with this dehnition. ^ 
The function U{x]k) is displayed in Fig. Elat t = 52, together with V{x;k) dehned by 
(j77|). for which the line of — estimated from the two-pion time correlator is also drawn. 
These hgures show that V (T; k) ~ —k^ and U (T; fc) ~ 0 for a: > i? within the statistical 
error for all quark masses. We hnd R ~ 10 (1.6 fm) for the heaviest quark mass, = 
0.736 GeV^. This stands for the largest R we obtained. This result signihes that the 
necessary condition for Liischer’s formula © is satished on the 24^ lattice for all our quark 
masses = 0.273 — 0.736 GeV^ with the current statistics of simulations. 


B. Scattering lengths 

We may estimate the scattering length from the wave function with two alternative 
methods : 

1. We extract k'^ by htting the asymptotic value of 17(x; k) to a constant and obtain the 
scattering length by substituting into (HD- The resulting /c^ and the scattering 
length are given in Table Q] in the column labeled with “from V”. We choose t = 52 
and the htting range Xm ^ x < -\/3L/2 (maximum value of x for lattice). The 
energy shift AVF is calculated from k"^ using ((221) • 

2. k‘^ is obtained by htting the wave function 0(T; k) with periodic Green function G{x\ k) 
dehned by taking and an overall constant as free parameters. An example of 
htting is illustrated in Fig. |21at t = 52 for = 0.273 GeV^, where the values from hts 
are shown with cross symbols. The htting range is the same as that for method^ The 
method for numerical evaluation of G{x] k) is discussed in Appendix El The ht works 
well, meaning that the contributions of Gim{x',k) and 6i{k) with I > 4 are negligible 
as expected. This point is conhrmed by htting with the function including G 4 o(T; k). 
The results are given in Table |l] (labeled “from 0”). 

We compare the scattering lengths obtained from the wave function with those from the 
conventional method of using the two-pion time correlator. We plot the normalized two-pion 
correlator R{t) of (j21(l in Fig. IHl which shows a clear signal that decreases exponentially in 
t. This means the absence of the singular-value solution for this volume, or otherwise R{t) 


^ The radius R thus defined has no direct relevance to the physical scale such as an effective range tq defined 
by k/ tan^o(fc) = l/oo + r^k^ -|- 0(fc‘^). It becomes larger as the statistical accuracy of and k) 
increases. We needed this somewhat artificial criterion to define R, unless otherwise we must appeal to 
some effective models. We also note that the rigorous estimation of the effects of the tail on the scattering 
length is not possible, unless Uk{x,y) for all x and y or the two-pion scattering amplitude off the mass 
shell for all energies is known. 
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should approach some constant. The effective masses for R{t) presented in Fig. [7| show the 
plateau over t ~ 22 — 70. Rows indicated with the label “from T” in Table Q] present the 
values of AW obtained by a single exponential htting for t = 24 — 68, /c^ estimated from 
AW, and the scattering length from Liischer’s formula fS) using thus estimated and 
goo{k) calculated by the numerical method discussed in ll|. We compare the scattering 
length obtained from three methods in Fig. |H1 The three methods give consistent results 
within statistical errors. We observe that the statistical errors for those from our new 
methods (“from V" and “from 0”) are signihcantly smaller than those from the two-pion 
time correlator (“from T”). This feature was experienced in the two-dimensional statistical 
models [l^. Our analysis so far is made at t = 52, but we checked that the results are 
independent of the choice of t in so far as t > 48. 

When one wants to obtain scattering lengths for the physical pion mass, there is yet an 
important problem of the chiral extrapolation. In Fig.|Hlwe carry out the chiral extrapolation 
of ao/rriTr using the data from method Q (he. “from R”) by assuming three htting forms : 


Fi{ml) 

F2{ml) 

Fsiml) 


2 


Ai -|- Bi ■ -\- Cl 

A 2 + B2/ml + C2-m; 

713(^1 + B 3 -ml- logiml/Cs] 


(30) 


The form of F 2 {rn^) is motivated from chiral symmetry breaking of the Wilson fermion and 
quenching effects, as suggested from quenched chiral perturbation theory [3,E3|- F^{rn^) 
is the form predicted by chiral perturbation theory (CHPT) for full QCD in the one loop 
order. Fi(m^) is a simple polynomial. We see that the three functions ht the data equally 
well and we cannot distinguish among them. The chiral limit of ao/m^(= Aj), however, 
depends sizably on the choice of the htting forms ; 


(1/GeV^) = —2.117(83) from Fi(m^) 
= —1.29(15) from F 2 {rri^) 
= —2.39(16) from F 2 ,{m‘l) 


(31) 


We cannot reduce this large systematic errors arising from the choice of htting forms, unless 
simulations are made close to the physical pion mass or the htting form is theoreti calh 
constrained. We add that the prediction from CHPT is an/mT, = —2.265(51) (1/GeV^) 


C. Results on small lattices 

We carry out the same analysis for 20^ and 16^ lattices to study the dependence on the 
hnite lattice size, with the numerical results presented in Tables im and ITm 

Our analysis on the 24^ lattice shows that the two-pion interaction range is at most 
R ~ 10 (1.6 fm) which happens for the heaviest quark mass, = 0.736 GeV^, so that 
Liischer’s formula is safely applied for the 24^ lattice for all our quark masses with the 
present statistical accuracy. This appears to indicate that 20^ is needed and 16^ may be too 
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small. This is not necessarily true, however, since R depends on the quark mass and the 
momentum which strongly depends on the lattice volume, as is seen by comparing the 
relevant entries in the three tables. 

To investigate the lattice size dependence of the interaction range, we plot V (T; k) and 
U{x; k) for the 20^ lattice in Fig. El and for the 16^ lattice in Fig. HOI both at t = 52. We 
cannot clearly observe a region x < L/2 where U{x]k) ~ 0 for heavier quarks, while such 
a region within statistics is visible for lighter quarks. The scattering lengths are calculated 
for the latter cases, i.e. at the two lighter quark masses = 0.273 and 0.351 GeV^ on the 
20^ lattice, and only at the lightest quark mass = 0.273 GeV^ on the 16^ lattice. 

Our compilation of the the scattering lengths, i.e. those obtained on the three lattice 
volumes for hve quark masses with three different methods, is depicted in Fig. Data 
points encircled by dotted lines are those from the two-pion time correlator for the case 
for which we cannot clearly hnd a region x < L/2 where the two-pion interaction vanishes. 
We do not hnd a signihcant volume dependence, however, for all quark masses including 
the case where the necessary condition for Liischer’s formula is not satished. The effects of 
deformation of the two-pion interaction due to hnite volume effects on AW of the two-pion 
system are apparently small compared with the statistical error. We emphasize, however, 
that the reliability of these data is guaranteed only after we obtain the results on the 24^ 
lattice where the necessary condition for Liischer’s formula is satished. 


V. CONCLUSIONS 

We have shown in this work that calculations of the two-pion wave function for the ground 
state of the 1 = 2 S'-wave two-pion system is feasible. We have investigated the validity of 
the necessary condition for Liischer’s formula and have found that it is satished for L > 24 
(3.92 fm) for = 0.273—0.736 GeV^. We have demonstrated that the scattering length can 
be extracted from the wave function with smaller statistical errors than from the two-pion 
time correlator which has been used in the studies to date. 

We have observed no signihcant volume dependence for the scattering lengths obtained 
from the two-pion time correlator at least for L > 16 (2.61 fm), in spite of the fact that 
the necessary condition for the formula is not satished in some cases for L = 16 and 20. 
The ehects of deformation of the two-pion interaction due to hnite size ehects on the energy 
eigenvalues of the two-pion are likely small compared with the statistical error. 

The present work opens the possibility to reduce statistical errors in the scattering phase 
shift with modest statistics. In this case, however, we have to concern with the contami¬ 
nation from inelastic scattering. This ehect is probably negligible for the ground state of 
the two-pion as in this work, but it may be important for the momentum excitation states, 
whose energies are close to the inelastic threshold. One may investigate ehects of inelas¬ 
tic scattering by evaluating I{x) in (jA2jl in Appendix fA1 from calculations of (j){x]k) and 
{'7i{'p)\7i{—x/2)\n7i;k). This would also clarify the error associated with our neglect of the 
inelastic channel, but the work is deferred to the future. 
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Another implication of the present work is the feasibility to calculate the decay width 
of the p meson through studies of the 1 = 1 two-pion system. While the evaluation of the 
disconnected diagrams with a good precision has been a computational problem, our new 
method, investigating the scattering system from the wave function of multi-hadron states 
in which the energy eigenvalue is extracted from the wave function at a single time slice, 
could lend a tactics that can be used to evaluate such complicated diagrams with a modest 
cost. 
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APPENDIX A: LUSCHER’S FORMULA FROM BETHE-SALPETER WAVE 


FUNCTION 


We discuss the derivation of Liischer’s formula m from Bethe-Salpeter (BS) wave func¬ 
tion in quantum field theory All considerations are made in Minkowski space, but these 
are not changed even in Euclidian space. 

We consider the BS function in the infinite volume defined by 


0oo(ufc) = (0|7ri(f/2)7r2(-f/2)|7ri(/c),7r2(-fc);in) , 


(Al) 


where we consider the particle as distinguishable and denote two distinguishable pions by tti 
and 7r2. The state |7ri(/c), 7r2(—/c); in) is an asymptotic two-pion state with the momentum 
k and —fc, and 7^i{x) and vr 2 (a;) are the interpolating operators for the pion tti and 7^2 at 
position X. Inserting complete intermediate energy eigenstates between the two fields, we 
obtain 

(j)oo{x]k) = J (-2^)32^ yZe^P'^A^7ri(^|7r2(-f/2)|7ri(fc),7r2(-fc);in) +/(f) , (A2) 

where Ep = I{x) is contribution of in-elastic scattering from the states with more 

than one pions, whose contributions are expected to be small for the energy 2Ek << 

This energy condition is satisfied for our case of the two-pion ground state, thus we neglect 
I{x) in the following. 

We decompose the BS function into the disconnected and the connected part, and rewrite 
them by the LSZ reduction formula as 


(7ri(^|7r2(-f/2)|7ri(fc), Ti 2 {-k); in) 


= (27r)^(5^(p — k) -I- 


H{p;k) 


2Epy/Ze^P'^/‘^ -k"^ -ie 

where H{p] k) is related to the off-shell two-pion scattering amplitude M(p; k) by 

H{p] k) = k) . 

^Pjptjk 

M{p] k) is defined from the pion 4-point Green function by 
M{p] k) 


(A3) 


(A4) 




—-1- — le 


- = d^z dV d^l/2 K(p,z)K(-ki,yi)K(-k2,y2) 


(0|T[7ri(z)7r2(x)7ri(yi)7r2(y2)]|0) , (A5) 


where all momenta and coordinates denoted by bold face characters refer to four-dimensional 
vectors and 


K(p,z) = 






(A6) 


In our case the momenta in (ESI) take 

ki = (Efc, fc) , k2 = (Efc, -k) 
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p = {Ep, p) , q = ki + ks - p = {2Ek - Ep, -p) . 


(A7) 


q is generally off-shell momentum and q = {Ep, —p) at on-shell {Ep = E^). The others are 
on-shell momenta. 

We assume that the scattering phase shift Si{k) with / > 1 is negligible, and regard 
H{p; k) and M(p; k) as functions of p = |p| and k = \k\. We also assume regularity for all p 
and k. M{p-, k) and H{p-, k) are normalized as 


at on-shell p = k. 
Substituting 


into 


(t)oo{x-,k) = 


M{k;k) = sin(5o(fc) 

k 

A 'jr 

H{k-,k) = —e*^“(^)sin5o(A;) 

K 


, we obtain 

d^p E[ (p; k) 


(AS) 


^ip-x 


ik 


(27r)3 p"^ — k‘^ — ie 


= e^^-^ + —H{t,k)Ukx) + P 

where we neglect an irrelevant overall constant. 

We assume 


d^p H{p;k)jo{px) 

{2'kY ~ 


(A9) 


h{x-, k) = 


d^p 


H{p; k)e-^P-^ 


. (27r)3- 

= —(A-|-fc) = 0 for a; > i? , 


(AlO) 


where h{x;k) depends only on x = |x| and k = \k\. In order to simplify (jA9jl . we use a 
formula : 

f d^p E{p) /■ ,, . k 


= J f{z)—no{kz) , 


(27r)3p2-P , ^ ^ , (All) 

where f{z) is an inverse Fourier transformation of E{p) which is a function of p = |p|. From 
this the third term in ()A9jl can be written by 

d^p H{p;k)jo{px) 


E{x-, fc) = P 


= J d^z g{z;k,x)^no{kz) , (A12) 


(27r)3 p2 — k'^ 

where g{z-, k,x) is an inverse Fourier transformation of i7(p; k)jo{px), which is given by 

d^p 


g(z\k,x) = 


(27r) 

dhlT 


-H{p;k)jo{px)e 


d^p 


H{p-, k)EP-^e-^P-^ 


J Att J (27r)^ 

= j ^ J 6 ^{y + x- z)h{y; k) , 

with the spherical coordinate O 3 . for x. By substituting HA 1. 'Ill into (lA 1 PH . we obtain 


(A13) 


E{x-,k) = J J d^y h{y,k)^no{k\x+ 


k 
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^ dy {ATiy‘^)h{y] k) J + 

POO U 

dy {^ny^)h{y;k) — 


Q{x-y) ■ no{kx)jo{ky) 
+e{y - x) ■ jo{kx)no{ky) 


(A14) 


where Q{x — y) is the step function, i.e. 0(x — y) = 1 for x — y > 0 and 0 for others. Under 
the assumption (lAlOll . (IA14|1 for x > i? can be written by 

k 

E{x;k) = —no(fcx) / dy {4:7iy‘^)h{y; k)jo{ky) 

4:71 Jo 


k 

= ^no(fcx) / dy {d'Ky^)h{y]k)jo{ky) 
47r Jo 

k 

= —no{kx)H{k]k) . 

47r 


(A16) 


Now we achieve the following simple expression of the BS function in the inhnite volume for 
X > R. 


ik 


k 


(j>oo{x', k) = + —H{k] k)jo{kx) + —H{k-, k)no{kx) 

= e 


47r 47r 

— io{kx) -\- cos5o{k)jo{kx) -\- sm.5o{k)no{kx) , (A16) 


where we use the relation (EHl) in the hnal step. 

The hrst term in ()A16f) is written in terms of ji{kx) as 


OO I 


Ak-x 


jo(kx) ^Y. {‘kx)i'i,{kx)Y7„(at)Yi„,{a,) , 


(A17) 


/=1 rn=—l 


with the spherical coordinate klk for k and Vt^ for x. Thus it does not contain the S-wave 
component and it is a regular function for all x. We hnd that the ratio of the coefficients of 
ioikx) and noikx) in (jAlbjl gives the S'-wave scattering phase shift. It should be noted that 
(f>oo{x; k) does not contain riiikx) with I > 1. This is attributed to the fact that we neglect 
the scattering phase shifts 6 i{k) with I > 1 and regard M(p; k) as a function of p = |p| and 
k = \k\. 

The BS function on a periodic box is dehned by 


(j){x-,k) = {0\7ii{x/2)7i2{-x/2)\7ii7i2;k) , (A18) 

where |7ri7r2; k) is the energy eigenstate of two-pion system on the periodic box with energy 
E = 2^771^ + /c^. Here we assume that the two-pion interaction range R is smaller than one 
half of the lattice extent, i.e. R < T/2, so that the boundary condition does not distort the 
shape of h(x; k). 

0(x; k) should be written in terms of 0oo(^; k) as follows. 

C50 I _ 

0(x; k) = Y.Y. k) , (A19) 

/=0 m=—l 
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where Cim are determined from the boundary condition together with the allowed value of 
k^, and k) is the Im component of 0oo(^; k) dehned by 

k)^ j k) . (A20) 

0(5;; k) satishes the the Helmholtz equation for x > R. The general solution of the 
equation on a periodic box Li^ can be written by (jSj) with Gim{x;k) dehned by ((7j). Thus 
the equation (iroi) yields 

OO I 

(p{x;k) = CimV^Yimikl^)(p^^ix;k) 

1=0 m=—l 

OO I 

= '^ VlmGlm{x]k) , (A21) 

;=0 m=—l 

for X > R. As mentioned in Sec. HH Gim{x]k) contains nii{kx) with only I' = 1 . (j)oc{Y, k) 
contains ni{kx) with only / = 0 as known from (lAlbfl and (IA17I1 . Thus only Uqo is non-zero 
and Gim{x] k) with / > 1 do not contribute in the second line of ()A21|1 . 

We use the expansion form of G{x] k) given by © to determine the allowed value of 
Gim and Vqq in (IA21(I . Comparing the S'-wave component of both lines of ()A21|1 . we hnd 

C'oo-e*''°('=)sin<5o(A;) = Uoo-^ (A22) 

dvr 

Coo ■ cos5o(^) = r’oo ■ 5'oo(^) • (A23) 


Finally we obtain Liischer’s formula m by taking the ratio of (IA22fl and ()A23fl . 


1 

tan 6 o{k) 


dvr 

T 


■ goo{k) . 


(A24) 


The other components of ()A2H) give only the informations for Gim- 

Let us make a comment on the condition (jAI Ojl . In quantum held theory this condition 
is generally not exactly satished and there can be a small tail in h{x-, k) for large x. Thus we 
cannot rigorously dehne the two-pion interaction range. Further an exact estimation of the 
ehects of the tail for the BS function is not possible, unless the oh-shell two-pion scattering 
amplitude H (p; k) for all p for given k is known. In this appendix we considered that the 
condition (jAIOjl is satished for some value of R, with a tacit assumption that the corrections 
from the interaction tails for the BS function are negligible. 


APPENDIX B: NUMERICAL CALCULATION OF PERIODIC GREEN FUNC¬ 
TION 


We rewrite G{x] k) in terms of dimensionless values as 


(dvr^L) ■ G(f; k) 


E 

m£'Z^ 


^im-X 
9 2 ’ 

— qf. 


(Bl) 
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where ql = /c^L^/(27r)^ Z) and X = x{27i)/L. The function G{x]k) can be expanded 
around the momentum = ni- (27r)^/L^ {rii G Z) as 


(dvr L) ■ G{x- k)= 2 _ „2 ^ 

P m^=nj 


2 _ „2 E o'"""' + T.(d- 4y-X(x-, j, n ,), 


i=i 


where 


F(X-,j,n^)= Y1 ■'e*™'-'' 


(B2) 


(B3) 


The function F{X]j,np) depends on the position X, lattice geometry and the expansion 
point Up. But it is independent of the physical quantity, such as quark masses and the 
strength of the two pion interaction. In our case of t he g round state of the two-pion, we set 
Up = 0. We can use the same techniques as in Ref. for the evaluation of the spherical 
zeta function. F{X] s, Up) takes hnite values for Re(s) > 3/2. The function at s = j G Z > 1 
is dehned by the analytic continuation from this region. 

First we divide the summation in F{X]s,np) into two parts as 


F(X-,s,n,)= E + E 


m?<inZ rrP>n^ 


(m^ - nl) * e* 


(B4) 


The second part can be written by an integral form as 


(m2e*™-^ 






m2>n2 ^ 

roo _ /*! _ rl _\ 

[ / '"F'e 

h “'0 2 “'0 ^jT[s 


—tim? —rip) im-X 


(B5) 


The hrst and second terms in ()B5jl converge at s = j G Z > 1, which are given by 




Afa-X 


E 


Am-X 


+ E7- 


J'- rrSF„2 ^“,2 ' /=! (j " ! „,^„2 " njY 


E 


g-(m2-n2) 


Am-X 


(B6) 


rrfi^ri^ '“P'' 


at s = j G Z > 1. The third term in (|B5l) is rewritten by Poisson’s summation formula ; 

J2nm-y (B7) 


E /(^) = E y d^y f{y)F 


and integration over y yields, 
t 


dt ^ 4_ Q-t{rn?-nl) ^im-X _ ^ ^-{X+2TTm)^ / (U) 


r(s) 


T(s)U/ 


^ e 

mgZ^ 


(B8) 


The hnal expression in (iBl converges at s = j G Z > 1 for X/{2x) ^ 7 ?. We do not need 
the values at X/(27r) = n G Z^, because these correspond to T = nL and these positions are 
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within the two-pion interaction range. Finally, gathering all terms and setting s = j 
in (B, we obtain 


F{X;j,np) = 


J 




E 




Am-X 


71 


3/2 


— r)! TT 5 [m^ 


nlY 


{j -1)\ Jo 


dt 


— (X+27rm)^/(4t) 


G Z > 1 


(B9) 
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FIGURES 



FIG. 1: Two-pion wave function (j){x;k) on 24^ lattice on {t,z) = (52,0) plane for = 
0.273 GeV^. The reference vector is set at xq = (7,5,2) (xq = |To| = 8.832). 
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FIG. 5; V{x;k) in unit of 1/a^ (left side) and U{x;k) (right side) on 24^ lattice at t = 52 for 
several quark masses. Horizontal axis is x = |x|. We plot a line at —k'^ obtained from the two-pion 
time correlator in left side of figures. 












































FIG. 6: Normalized two-pion correlator R{t) on 24^ lattice for = 0.273 (lightest), 0.44 and 
0.736 GeV^ (heaviest). Scale of vertical axis is log-scale. 
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FIG. 7: Effective mass for R{t) on 24^ lattice in unit of 1/a at ml = 0.273 (lightest), 0.444 and 
0.736 GeV^ (heaviest). 
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m/(GeV2) 


FIG. 8; Pion scattering length ao/m 7 r (1/GeV^) on 24^ lattice obtained from the two-pion time 
correlator (“from T”), the constant fit of V{x]k) (“from V") and the fitting the wave function 
(“from (/>”). Results of the chiral extrapolation of ao/m-Tr obtained from V(x;k) (“from F”) with 
the three fitting functions are also plotted lines. 
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FIG. 9; V{x;k) in unit of 1/a^ (left side) and U{x;k) (right side) on 20^ lattice at t = 52 for 
several quark masses. In left side of figures we also plot a line at —obtained from the two-pion 
time correlator. Horizontal axis is x = Ixl. 
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FIG. 11; Volume dependence of pion scattering length (1/GeV^) for five quark masses 

obtained from the time correlator (“from T”) and the wave function (“from V” and “from (/>”). 
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TABLES 


L 3 = 243 


mjr (GeV) 

0.85773(27) 0.76675(28) 0.66649(30) 0.59205(31) 0.52290(33) 

ml (GeV^) 

0.73570(47) 0.58790(43) 0.44421(40) 0.35052(37) 0.27342(35) 

Xm 

9.8 

9.8 

9.8 

9.4 

9.0 

AVL (xlO-3 GeV) 






from T 

1.748(64) 

2.011(73) 

2.337(86) 

2.60(10) 

2.85(11) 

from V 

1.817(25) 

2.049(29) 

2.336(35) 

2.559(37) 

2.804(40) 

from (j) 

1.832(43) 

2.071(49) 

2.346(58) 

2.550(65) 

2.784(66) 

k^ (xlO-3 GeV^) 






from T 

1.500(55) 

1.543(56) 

1.559(58) 

1.540(59) 

1.492(60) 

from V 

1.559(22) 

1.572(22) 

1.558(23) 

1.517(22) 

1.468(21) 

from (j) 

1.572(36) 

1.589(38) 

1.565(39) 

1.512(39) 

1.457(35) 

ao/m^ (1/GeV^) 






from T 

-0.969(32) 

-1.112(36) 

-1.291(42) 

-1.439(49) 

-1.585(57) 

from V 

-1.003(12) 

-1.131(14) 

-1.291(17) 

-1.420(18) 

-1.562(20) 

from (j) 

-1.010(21) 

-1.141(24) 

-1.296(29) 

-1.415(32) 

-1.552(33) 


TABLE I: Results on 24^ lattice from two-pion time correlator (“from T”), constant fit of R(x; k) 
(“from R”) and fitting wave function with G{x;k) (“from (/>”). The fitting range of V{x]k) and 
4>{x;k) are quoted as Xm- 
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= 20^ 


(GeV) 

0.85754(25) 0.76653(26) 0.66628(28) 0.59188(30) 0.52283(32) 

(GeV^) 

0.73537(42) 0.58757(40) 0.44393(38) 0.35033(36) 0.27335(34) 

Xm 

9.0 8.2 

AIT (xlO-3 GeV) 

from T 

3.219(88) 3.66(10) 4.20(11) 4.63(13) 5.03(15) 

from V 

4.457(63) 4.928(68) 

from (j) 

4.62(12) 5.05(12) 

k^ (xlO-3 GeV^) 

from T 

2.763(76) 2.810(77) 2.805(76) 2.744(75) 2.635(78) 

from V 

2.643(38) 2.583(36) 

from (p 

2.741(72) 2.646(64) 

ao/rriT, (1/GeV^) 

from T 

-1.003(24) -1.139(27) -1.310(31) -1.447(35) -1.584(41) 

from V 

-1.401(17) -1.556(19) 

from (p 

-1.446(33) -1.589(34) 


TABLE II: Results on 20^ lattice from two-pion time correlator (“from T”), constant fit of V(x; k) 
(“from E”) and fitting wave function with G{x;k) (“from (p'”). The fitting range of V{x]k) and 
(p{x; k) are quoted as Xm- 
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= 103 


(GeV) 

0.85759(34) 0.76672(37) 0.66660(41) 0.59225(45) 0.52318(50) 

(GeV^) 

0.73546(58) 0.58786(57) 0.44435(55) 0.35076(53) 0.27371(52) 

Xm 

7.8 

AIT (xlO-3 GeV) 

from T 

6.51(17) 7.28(20) 8.20(23) 8.93(26) 9.62(32) 

from V 

9.64(21) 

from (j) 

9.71(41) 

k^ (xlO-3 GeV^) 

from T 

5.59(15) 5.59(15) 5.48(15) 5.31(16) 5.06(17) 

from V 

5.07(11) 

from (j) 

5.10(22) 

ao/m-„ (1/GeV^) 

from T 

-1.002(22) -1.122(26) -1.270(30) -1.393(35) -1.517(43) 

from V 

-1.519(28) 

from (j) 

-1.528(55) 


TABLE III: Results on 16^ lattice from two-pion time correlator (“from T”), constant fit of V(x; k) 
(“from V”) and fitting wave function with G{x;k) (“from cp”). The fitting range of V{x]k) and 
(p{x; k) are quoted as Xm- 
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